Astronomy & Astrophysics manuscript no. articolo-processo-s 


February 2, 2008 


(DOI: will be inserted by hand later) 





The s-process weak component: uncertainties due to convective 

overshooting 

V. Costa 1 , M. L. Pumo 2 , A. Bonanno 1 , and R. A. Zappala 2 

' Istituto Nazionale di Astrofisica - Osservatorio Astrofisico di Catania (INAF-OAC), Citta Universitaria, Via S. Sofia 78, 
1-95123 Catania, Italy 

2 Universita di Catania, Dipartimento di Fisica e Astronomia (sez. astrofisica), Citta Universitaria, Via S. Sofia 78, 1-95123 
Catania, Italy 

Received 1 August, 2005; accepted 29 September, 2005 

Abstract. Using a new s-nucleosynthesis code, coupled with the stellar evolution code Star2003, we performed simulations 
to study the impact of the convection treatment on the s-process during core He-burning in a 25 M Q star (ZAMS mass) with 
an initial metallicity of Z=0.02. Particular attention was devoted to the impact of the extent of overshooting on the s-process 
efficiency. The results show enhancements of about a factor 2-3 in s-process efficiency (measured as the average overproduction 
factor of the 6 s-only nuclear species with 60 < A < 90) with overshooting parameter values in the range 0.01-0.035, compared 
to results obtained with the same model but without overshooting. The impact of these results on the p-process model based on 
type II supernovae is discussed. 
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1. Introduction 

As already pointed out by Clayto n~(1968)| it is now widely ac- 
cepted that s-nuclei are formed via neutron exposures on iron- 
peak nuclei, and it is also clear that more than one event (or 
s-process "component") with a single set of physical condi- 
tions (like neutron exposure t, initial abundances, and neutron 
density) is necessary in order to account for the observed solar 
distribution of s-nuclei abundances. 

Current views on the subject suggest the existence of two 
components. One is the so-called "main" component that is as- 
sociated with low-intermediate mass stars (~ 1-5 M Q ) during 
their AGB phase and that should be responsible for the synthe- 
sis of s-species from A - 90 up to 209 Bi (Gallino et al. 119981 
Busso et aLEHHEOnil Herwig et al. l2U0H Lugaro et al. l2U0H 
Goriely and Siess 2004 1. The other is the so-called "weak" 
component that covers the 60 < A < 90 mass range, and is 
found in massive (M > 10 M Q ) stars during their core He- 
burning phase, as first proposed by Peters (1968 1. Moreover, 
a particular class of massive stars that may contribute to the 
s-process enrichment of interstellar medium includes the so- 
called Wolf-Rayet stars (Meynet et al. 120011 and references 
therein). Thanks to the loss of their initial H envelopes, these 
stars may show the products of He-burning on their surface. 

Some authors (Gallino et al. ITWgl Busso et al. [T333] 
Kappeler fTgggi Goriely and Siess |2"0UT1 12004 Lugaro et al. 
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2003) also suggest the existence of a "strong" component that 
should be responsible for the synthesis of s-species around 
208 Pb. This component should occur in low-metallicity stars of 
low-intermediate mass during the AGB phase. Observational 
evidence of an enhanced lead abundance in metal-poor stars 
(Van Eck et al. 2001 ) seems to confirm this possibility. 

Although the general features of s-nucleosynthesis seem to 
be well-established, there are still ambiguities and uncertainties 
concerning the weak component. Some authors (see e.g. Raiteri 
et al.rE?9H The et al.l2"U0ll Hoffman et al. 1277011 Woosley et 
al. 2002 1 suggest that this process might be prolonged during 
subsequent post-He-burning stellar evolutionary phases, so that 
the s-abundances at the end of core He-burning could be mod- 
ified, but this hypothesis needs further investigation (Arcoragi 
et al. fTWTl The et al. l2UUUl Rayet & Hashimoto 2000 Woosley 
et al. 2002). Besides, both the nuclear mass range listed above 
for the weak component and the overall efficiency of the pro- 
cess are subject to significant uncertainties that come from both 
the nuclear physics and astrophysical modelling aspects of the 
problem. 

Uncertainties due to nuclear physics are connected to the 
value of the l2 C(a, y) 16 reaction rate, with its impact on the 
stellar core structure, and to reaction rates on which the so- 
called "neutron economy" (the balance between neutron emis- 
sion and captures) is based. Many authors (see e.g. Rayet & 
HashimotoEOnni The et al. l^UOT)! Rayet et al. l2TKJTl Costa et al. 
2000 2003 hereafter Paper I and Paper II) have examined the 
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impact of these uncertainties on the s-process weak-component 
yields. On the astrophysical side, the lack of a self-consistent 
convection theory is still one of the major problems in calculat- 
ing stellar structures. Thus it can be a significant source of un- 
certainty for nucleosynthesis processes, like the s-process weak 
component. 

Some weak component s-process calculations performed 
so far have been developed using stellar evolution codes that 
are based on simple mixing length theory without overshooting 
(see e.g. Rayet & Hashimoto 2000, The et al. 2000). However, 
overshooting may affect s-process nucleosynthesis in massive 
stars by giving rise to an increase in the convective core mass. 
This may affect the s-process through lengthening of the He- 
burning lifetime and an increase in the amount of material that 
experiences neutron irradiation. Indeed, matter originally lying 
in stellar regions with temperatures that are too low to allow the 
main neutron source to work ( 22 Ne(a, n) 25 Mg becomes effec- 
tive at temperatures T > 2.5 x 10 8 K) would be mixed towards 
central hotter regions, while a change of the He-burning life- 
time would directly change the s-process lifetime. 

Langer et al. (1989) used two 30 M Q star models (with- 
out and with overshooting) to study the overshooting effect on 
the s-process in massive stars. They conclude that overshoot- 
ing leads to an increase in s-nuclei production (see their Table 
2) with an enhancement factor of about 1.5-2 on final mass 
fractions and an increase in the average number of captured 
neutrons per nucleus of 56 Fe from a value of 3.6 to a value of 
4.5. 

We re-examined the impact of overshooting on the s- 
process not only because a comprehensive study of the role 
of overshooting on s-nucleosynthesis is needed in the frame- 
work of recent nuclear data and stellar models, but also because 
the approach of Langer et al. ( 1989 1 for overshooting is based 
on a simple artificial enhancement of the convective core ex- 
tension by a 40% fraction of pressure scale height, while we 
used a more sophisticated approach based on a diffusive algo- 
rithm. In this approach an "overshooting parameter" / deter- 
mines the extension of convectively mixed core, with the over- 
all efficiency of overshooting described by a diffusion equa- 
tion (see Sect. l2.1l for more details). This approach has already 
been tested with some success for the Sun in order to find a 
solution to the lithium depletion puzzle (see e.g. Blocker et al. 
119981 Schlattl & Weiss I1999> . Of course the overshooting in 
massive stars is different, as it is concerned with the "outer" 
border of the central convection zone, while the Sun is charac- 
terised by overshooting through the "inner" border of the ex- 
ternal convection zone. However, Salasnich et al. ( 1999 1 inves- 
tigated the internal convection zone in massive stars using dif- 
fusion to model the convective overshooting. They give some 
indications of how to modify the overshooting parameter in or- 
der to account for the observed distribution of massive stars 
across the HR diagram, and, in particular, they claim a value 
of about / = 0.015 for 20 M Q stellar models for a setting of 
this parameter. Other authors (Young et al. 2001 1 suggest that 
the overshooting parameter may be mass-dependent and that 
massive stars probably require a higher / value. Thus, since 
the proper / value is still being debated, we examine a whole 
range of its possible values in this paper. 



Moreover, our study may be of interest in p-process mod- 
ellings, because some of the problems in the model for the p- 
process (taking place in the type II supernovae O-Ne layers) 
might be linked to uncertainties in the s-process weak compo- 
nent efficiency (Paper I; Paper II; Rayet et al. 2001 Arnould & 
Goriely 2003 1, since the relevant nuclides are p-process seeds. 

2. Input physics 

2.1. Stellar model 

A preliminary set of nucleosynthesis calculations was per- 
formed with the same stellar data as in Papers I and II. It is 
basically the model of a He star of 8 M Q taken from Nomoto 
& Hashimoto (1988 1, which is supposed to correspond to a 
25 M at the ZAMS. This model treats convection with the 
standard mixing length theory and does not include overshoot- 
ing or mass loss. 

Our new calculations were performed with the stellar evo- 
lution code Star2003, with the Garching stellar evolution code 
(Weiss & Schlattl 2000 1 as its basis, and the models were 
evolved starting from ZAMS until the end of core He-burning. 

Mixing was modelled within a diffusion approach (Blocker 
et al. 1 19981 Freytag et al. 1996), and abundance changes of the 
12 considered nuclear species were calculated with the follow- 
ing diffusion equation: 



dX _ (dX\ d 

dt \ dt Inuc dm <- 



K") ~dm r 



(1) 



where the first term on the right is the time derivative of a given 
nuclidic mass fraction due to nuclear reactions, while the sec- 
ond is the diffusive term that describes mixing. The derivatives 
were calculated with respect to the radial stellar mass m r , while 
the other symbols used in Eq. [Qhave their standard meaning. 
The difference among convective, overshooting, and radiative 
regions lies in the value used for the diffusion coefficient D. In 
convective zones, the diffusion coefficient is given by: 



D, 



(2) 



where v c is the average velocity of convective elements derived 
from the mixing length theory, and I is the mixing length. For 
the overshooting regions, the following diffusion coefficient 
was used instead: 



D omr = D exp-^; 

tly 



H v =f- H p 



(3) 



where Do is taken as equal to the value of D com , at the up- 
per radial edge of the convection zone established through the 
Schwarzchild criterion; z = \r - r ec i ge \ is the radial distance from 
the same edge; H p is the pressure scale height, while / is the 
overshooting parameter that determines the overall efficiency 
of overshooting and the extension of the convectively mixed 
core. For radiative zones, the diffusion coefficient was put to 
zero, and abundance changes were only due to the nuclear re- 
action term. 

As for energy generation and chemical evolution, the code 
includes a hydrogen burning network (pp chain, CNO cycle) 
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and a He-C-O network, following the evolution of 12 nuclides 
CU, 3 He, 4 He, 12 C, 13 C, 14 N, 15 N, I6 0, I7 0, 20 Ne, 24 Mg, 28 Si). 
Reaction rates were taken from Caughlan et al. (Q985l. The 
rate of l2 C(a, y) 16 is about 20% higher than the one reported 
in the NACRE compilation (Angulo et al. 1999 1, which is also 
the one included in the s-process reaction network described 
in Sect. 12.2.11 This choice was made in order to make our 
results comparable to those obtained in Papers I and II. Note 
that the 12 C(a, y) 16 rate from Caughlan et al. ( 119851 is inside 
the uncertainty range reported in NACRE. Calculations entirely 
based on updated rates will be presented in a future paper. 

As for the equation of state, the OPAL01 EOS tables are 
used, which are an updated version of the OPAL96 tables, 



revised according to the prescriptions by Rogers (2001) (see 
Bonanno et al. 2002 for details). The OPAL opacity (Iglesias 
& Rogers 1996) are used, implemented in the low-temperature 
regime by the tables of Alexander & Ferguson ( 1994). 



2.2. The nucleosynthesis code 

Our new s-nucleosynthesis code uses a differential equation 
solver that is based on the implicit Kaps-Rentrop method, 
which provides an estimate of the "truncation error" associ- 
ated with each single integration step (see Press et al. 1996 for 
details). In this way settings concerning precision requests can 
be based on truncation error constraints. Some numerical tests 
on the Kaps-Rentrop method and a comparison with other nu- 
merical solving procedures can be found in Timmes (1999). 
The s-process calculations reported by Papers I and II, Rayet 
et al. {1995}, and Rayet & Hashimoto (2000) were based on 
the Wagoner two-step linearisation procedure (Wagoner ll969> . 
The Wagoner method does not provide any error estimate, so 
that precision is only based on putting a limit on the relative 
variations of isotopic abundances from one time step to the 
next. 

The "post-processing" technique is used to couple the nu- 
cleosynthesis simulations with stellar evolution data. This tech- 
nique can easily be applied to radiative zones, where nucle- 
osynthesis occurs according to local physical conditions. The 
prescriptions described by Prantzos et al. ( 1987) are used for 
the convective core, but as they do not say how to include the 
effects from changes in the convection zone borders' position 
with time, we added the prescription described by Kippenhahn 
& Weigert (I1990> to include this effect. 

If X is the "local" rate of change of the mass fraction of 
a given nuclide (due to local physical conditions), ;«i and 1112 
are the time derivatives of convection zone borders' position 
(mass coordinates), M c = (ni2 - mi) is the convection zone 
mass extension, and X\ and X2 the mass fractions of the same 
nuclide just below and above the convection zone, as a mass 
fraction discontinuity can generally be present in convection 
zone borders. Then if^- is the average mass fraction uniformly 
valid over the whole convection zone for the given nuclide, then 
the rate of change of x may be expressed as follows: 



The first term on the right needs the evaluation of "averaged" 
rates considered by Prantzos et al. (1987), while the second 
and third terms were implemented in our code in order to treat 
changes in the convection zone borders. The second term on 
the right side of Eq.|4]is considered only when n\\ < 0, while 
the third term is considered only when n%2 > 0, which means 
that these two terms take into account the inclusion of material 
that is not previously mixed into the convection zone when it 
expands. That is why the two terms above are not considered 
when the convection zone shrinks. 

Another aspect that has to be taken into account when cou- 
pling the stellar model with the s-nucleosynthesis code is re- 
lated to the mixing inside the convective core. While the stel- 
lar evolution code uses diffusion to model the mixing of the 
species involved in the 12-nuclide network described above 
(Sect. 12. 1> . instantaneous mixing is assumed for all s-nuclides 
where a local equilibrium is assumed to hold for the neutrons 
(see Prantzos et al. 1987 for details). Thus the overall procedure 
is not strictly self-consistent; but we believe that, while the in- 
clusion of diffusion in the stellar code allows us to model the 
stellar core structure evolution better (in terms of temperature 
profiles, mass extension, He-burning lifetime), neglecting dif- 
fusion in the nucleosynthesis code should not have important 
consequences, as long as the wide uncertainties in stellar struc- 
ture due to convection uncertainties are under study. A more 
refined and self-consistent approach will be presented in a fu- 
ture paper. 

2.2.1. The reaction network 

The s-nucleosynthesis code uses a much wider nuclear net- 
work than the one used by the stellar evolution code. This 
network is the same as the one in Papers I and II, Rayet & 
Hashimoto J2000t . and Rayet et al. (200T). It includes 472 
nuclides and 834 reactions. Differences are present in the re- 
action rates; in particular, rates from the Network Generator 
NETGEN 1 (Aikawa et al. l2005l were used. 

In order to test the dependability of the s-nucleosynthesis 
code, a preliminary s-process calculation was performed with 
the same rates and stellar model as in Papers I and II, since our 
results (see Table reproduced theirs quite well; we refer to 
the results obtained by them with NACRE adopted values for 
22 Ne(o', n) 25 Mg and 22 Ne(or, y) 26 Mg reaction rates. 



3. Models 

To estimate the impact of convective overshooting on the 
s-process efficiency, we calculated the evolution of a Z = 
0.02, 25 M Q star 2 with the code Star2003 and with mix- 
ing length parameters a = 1.7 and 2.0 and, for overshoot- 
ing parameters, / = 10~ 5 (model without overshooting), 
0.01, 0.015, 0.02, 0.025, 0.03, and 0.035. 

The s-process efficiency can be analysed in terms of the 
following parameters: 



X = 



1 



f 

mi 



Xdm + (x - X\)rii\ + (X 2 - x)mi 



(4) 



1 Available online at http://www.astro.ulb.ac.be/Netgen/ 

2 Initial composition was taken from Anders & Grevesse ( 1989 1 
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the average overproduction factor Fq for the 6 s-only nuclei 
70 Ge, 76 Se, 80 Kr, 82 Kr, 86 Sr, and 87 Sr given by 



with Ft 



where F , is the overproduction factor, X,- the mass fraction 
(averaged over the convective He-burning core) of s-only 
nucleus i at the end of s-process, X, ,„, the initial mass frac- 
tion of the same nucleus, and N s — 6 is the number of s-only 
species considered; 

the average overproduction factor Fq identical to Fq, but for 
the 18 s-nuclei 3 within the mass range 88 < A < 130; 
the maximum mass number A max , for which the species in 
the 60 < A < A max mass range are overproduced by at least 
a factor of about 10 and 5 (respectively, the first and second 
value in the fourth column of Tables ^and[2] and first and 
second value in the third column of Table 0; 
the number of neutrons captured per 56 Fe seed nucleus n c , 
defined as 

: "" [Y A - K ,(())| 



= £ (A -56)1 



A=57 



*56(0) 



where Y A is the mole fraction of the species with mass num- 
ber A at the end of the s-process, Y A (Q) the initial mole frac- 
tion of the same nucleus, and ^(O) the initial 56 Fe mole 
fraction (the abundances are averaged over the convective 
He-burning core); 

- the maximum convection zone mass extension (hereafter 
m.c.z.m.e.) during the core He-burning s-process; 

- the duration of core He-burning s-process (value in the 
sixth column of Tables ^ and |2| and value in the fifth col- 
umn of Tabled. 

The results of our computations are summarised in terms of the 
previous parameters in Tables ^ an d ID f° r the mixing length 
parameters a — 1.7 and 2.0, respectively. The results of s- 
process test simulation performed with data used in Papers I 
and II are reported in Table [3] Figures ^ and|2| show the over- 
production factors of the 6 s-only nuclei within the mass range 
70 < A < 90 for a = 1.7 and 2.0. 

4. Results 

The data reported in Tables ^ an d 121 show the following main 
features: 

1. a sudden enhancement of the s-process efficiency when 
passing from a "no-overshooting" model (/ = 10~ 5 ) to 
any model that includes overshooting, with the following 
enhancements for the s-process indicators: and Fq en- 
hanced by a factor ~ 2-3 and ~ 1.5-2, respectively; n c in- 
creased from values of 4.2-4.4 to values of 5.2-5.6; A max 
enhanced from a value around 90 to values of 104-1 10, at 
least if a minimum overproduction factor of 5 is used to 
define A max . Even the minimum value adopted for the over- 
shooting parameter (/ = 0.01) is enough to produce a sig- 
nificant increase in all the adopted s-process indicators. 



0.035. 

0.03..; 

0.02. 

0.01.. 
0.025.. 
0.015..- 




Fig. 1. Overproduction factors of the 6 s-only nuclei within the 
mass range 70 < A < 90 as a function of mass number for stel- 
lar models with mixing length parameter a = 1.7. The labels 
refer to the overshooting parameter value. 




75 80 
mass number 



3 These are the nuclei indicated as "mainly" produced by the s- 
process in Anders & Grevesse 1 1989 1. 



Fig. 2. Same as Fig.^but for stellar models with mixing length 
parameter a = 2.0. 



a slow change of s-process efficiency when passing from 
one / value to the next. This is shown more clearly in Figs. 
\l\ and |3 where the overproduction factors for the lightest 
six s-nuclei are shown for the different / values. For in- 
stance, the biggest change in Fq is about a 50% enhance- 
ment from the lowest to the highest value for each of the 
two data sets (a = 1.7 and a = 2.0). 

some anomalies concerning the pattern of the s-process effi- 
ciency, which does not seem to grow monotonically strictly 
with the chosen / values. For each of the two data sets, the 
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Table 1. Values of the "indicators" of the s-process efficiency (see text) for stellar models with mixing length parameter a - 1.7. 
The overshooting parameter value used in the stellar evolution code is reported in the first column 



/ 


Fo 


^0 


4 

' 1 max 


n c 


m.c.z.m.e. 


Duration [sec] 


io- 5 


99.88 


6.57 


91 - 


96 


4.24 


5.31M 


2.43 


•10" 


0.01 


246.13 


11.16 


94- 


104 


5.22 


6.39M e 


2.27 ■ 


10 13 


0.015 


216.74 


10.13 


94- 


104 


5.09 


6.87M 


2.09' 


10 13 


0.02 


257.81 


11.57 


94- 


104 


5.32 


7.23M S 


2.06 


10 13 


0.025 


236.40 


10.82 


94- 


104 


5.20 


7.31M Q 


2.16 ■ 


10 13 


0.03 


304.17 


13.20 


96- 


104 


5.55 


7.97M G 


2.48 ■ 


10 13 


0.035 


310.39 


13.49 


96- 


104 


5.58 


8.67M 


2.44 ■ 


10 13 



Table 2. Same as Table^but for stellar models with mixing length parameter a = 2.0. 



/ 


^0 


F 


4 


n c 


m.c.z.m.e. 


Duration [sec] 


io- 5 


136.74 


7.41 


91 - 


96 


4.48 


5.33M 


2.40 ■ 10 13 


0.01 


331.44 


14.29 


96- 


110 


5.61 


6.59M 


2.32 • 10 13 


0.015 


222.21 


10.28 


94- 


104 


5.13 


6.74M G 


2.11 • 10 13 


0.02 


259.42 


11.57 


94- 


110 


5.31 


7.16M G 


2.07 • 10 13 


0.025 


264.23 


11.81 


94- 


110 


5.36 


7.38M e 


2.09 • 10 13 


0.03 


309.26 


13.35 


96- 


110 


5.59 


7.93M e 


2.13 • 10 13 


0.035 


316.68 


13.68 


96- 


110 


5.64 


8.63M G 


1.97 • 10 13 



Table 3. Same as Table^but for calculations performed with the data used in Papers I and II. 



Fo 


Fo 


4 


n c m.c.z.m.e. 


Duration [sec] 


76.89 


4.06 


89-91 


3.84 5.66M Q 


2.42 ■ 10 1J 



main anomaly involves the model with / = 0.01, which 
gives birth to a more efficient s-process than the one ob- 
tained with higher / values. A minor anomaly is noticed 
for the model with / = 0.02 and a — 1.7 as well. 

5. Discussion 

5.1. Impact of overshooting 

The above results clearly show the level of uncertainty in the 
modelling of the weak s-process component due to the cur- 
rent lack of a self-consistent theory of stellar convection. More 
specifically, the current uncertainties due to overshooting actu- 
ally do give rise to a factor ~ 2-3 uncertainty in the s-process 
efficiency, an uncertainty that is linked more to the presence 
or absence (/ = 10~ 5 ) of overshooting than to the exact value 
of the overshooting parameter. This situation would perhaps 
change in post-He-burning phases. 

The not strictly monotonic increase in the s-process effi- 
ciency with the / value (point 3. Sect. @}, may be due to the 
non-linear nature of the stellar structure and evolution equa- 
tions. Looking at the values reported in Tablesn an d|2] it ap- 
pears that the m.c.z.m.e. (and the convection zone mass exten- 
sion, in general) grows monotonically with /, but the same is 
not true for the He-burning phase lifetime. It can be noticed that 
the pattern followed by the He-burning lifetime values roughly 
follows the one of the Fq and n c values for both data sets. This 
could be interpreted as follows: the extension of the convection 
zone has a direct impact on the amount of material available 
for He-burning and consequently on the s-process duration and 
neutron exposure. Thus, the longer duration of the s-process 



for the models with / = 0.01 (compared to other models with 
overshooting) is probably the reason for the enhancement of 
s-process efficiency. 

Other anomalies are present in our data. The efficiency of 
the a — 1.7, / = 0.025 model is lower compared to the one 
of the a = 1.1, f = 0.02 model, despite the fact that corre- 
sponding lifetimes follow a "correctly" increasing order. This 
may be due to a peculiarity of this stellar model (see Fig. [3Ji: 
the / = 0.025 model has a higher m.c.z.m.e., but for most of 
the time the / = 0.02 model shows a larger convection zone. 
The s-process efficiency reflects this feature, too. 

5.2. Possible contribution from subsequent burning 
stages 

Arcoragi et al. (1991 1 suggest that core C-burning and shell 
He-burning should not cause significant changes to the heavy 
(A > 60) nuclide abundances, but the models they developed 
do not include overshooting. Other authors (Raiteri et al. 119931 
The et al. l2000l Hoffman et al. l200T1 Woosley et al. l2002l sug- 
gest that shell C-burning and other subsequent phases might 
significantly modify the heavy nuclide abundances. According 
to Tables^andlJKsee the m.c.z.m.e. value), the adopted / value 
does have a big impact on the maximum extension of the con- 
vection zone, and a value near 9 M is obtained for / = 0.035. 
This might bring s-nuclei that are already produced and some 
22 Ne towards layers that could be engulfed into subsequent 
shell-burning stages. Some additional investigation of the pos- 
sible contribution from post-core-He-burning stages should be 
performed, taking uncertainties from stellar evolution and nu- 



6 



V. Costa, M. L. Pumo, A. Bonanno, and R. A. Zappala: The s-process weak component 




5.0X10 12 1.0x1c 15 1.5x10 13 2.0X10 13 2.5x10 13 



time [sec] 

Fig. 3. Mass extension of the convection zone (in unit of M ) as 
a function of time during the core He-burning s-process for the 
a = 1.7, / = 0.02 model (solid line) and for the a = 1.7, / = 
0.025 model (dotted line). 

1 oooo p ' 1 1 1 1 ' ' 1 ' 1 1 1 ' 1 1 = 




60 80 100 120 140 

Fig. 4. Overproduction factors for the s-nuclei within the mass 
range 70 < A < 130. Solid lines show results coming from 
Paper I, while the dotted line shows our results for the a = 
2.0, / = 0.01 model. The labels fl, (i=l to 3) refer to the value 
of the 22 Ne(a, n) 25 Mg reaction rate used for s-process calcula- 
tions in Paper I (see Paper I for the definition of the various R, 
values). 

clear physics into account. We are convinced that uncertain- 
ties in stellar evolution calculations, particularly those concern- 
ing post-He burning phases, do not allow any final conclusion 
about the possible contribution of post-He burning phases to 
the s-process yields. 



5.3. Impact on the p process 

A comparison of our s-process results with those of Paper I 
(see Fig.|4j shows that our most efficient s-process is situated 
between the results obtained in Paper I with the rate Ri and the 
ones with the rate R2, which means that our results only seem 
to help a little for the well-known problem of the underproduc- 
tion of light p-isotopes in current p-process models (Paper I, 
Rayet et al. 2001 1. However, given the fact that the argument 
for the possible enhancement of the 22 Ne(a, n) 25 Mg reaction 
rate is still being discussed (Jaeger et al. 2001 Koehler 2002 1, 
a combination of a "moderately" enhanced s-process (referring 
to what would be needed for the p-process model to be satis- 
factory) by means of a stellar model based on overshooting, 
and a slightly enhanced 22 Ne(o', n) 25 Mg reaction rate (like the 
R2 value reported in Paper I) might possibly bring the light p- 
isotopes abundances coming from the p-process type II super- 
nova model in line with the solar distribution. 



6. Summary and final comments 

Many s-process simulations have been performed so far by var- 
ious authors, but there is still a lack of a detailed scrutiny of 
the impact of stellar evolution modelling uncertainties on the 
s-process weak component. In this paper we have shown the 
role of overshooting on the s-process in massive stars. 

Stellar models of 25 M with initial Z = 0.02 metallicity 
have been developed until He exhaustion in the core, using a 
range of values for the overshooting parameter /. Then models 
were used to "post-process" the stellar nuclidic composition 
with a new s-nucleosynthesis code. 

Models with overshooting give higher s-process efficiency 
compared to no-overshooting models, as concluded by Langer 
et al. (1989 1 with their simplified treatment of overshooting. 
However, our enhancements in the average overproduction 
(factors of about 2 to 3) are slightly higher than theirs (fac- 
tors of about 1 .5 to 2), and similar behaviour is found for the 
corresponding n c values (4.2-5.6 and 3.6-4.5, respectively), al- 
though their stellar model treats a 30 M star, which should be 
expected to give a more efficient s-process than our 25 M stel- 
lar models (Prantzos et al. [19901 The et al. l2000T l. 

Less important, but not negligible, variations are obtained 
when changing the overshooting parameter / value in the 
0.01 - 0.035 range, and the link between the / value and the 
s-process indicators values is nearly monotonic. 

As regards the connection between s- and p-process, the 
increase in efficiency of the s-process might possibly help with 
the problem of underproduction with respect to solar of light 
p-isotopes in type II supernovae. 

Other studies of stars with different mass and metallicity 
are needed to confirm our conclusions. 
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